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WroNG - Wroclaw Neutrino Generator of events for single pion 
production 

Jan T. Sobczyk*^* , Jaroslaw A. Nowak^, Krzysztof M. Graczyk'* 

'^Institute of Theoretical Physics, Wroclaw University, 
pi. M. Borna 9, 50-204 Wroclaw, Poland 

We constructed a new Monte Carlo generator of events for neutrino CC single pion production on free nucleon 
targets. The code uses dynamical models of the DIS with the PDFs modified according to the recent JLab data 
and of the A excitation. A comparison with experimental data was done in three channels for the total cross 
sections and for the distributions of events in invariant hadronic mass. 



1. INTRODUCTION 



Neutrino-nucleon (or nucleus) CC interactions 
are described in a framework of three different 
theoretical schemes: quasi-elastic, resonance, and 
deeply inelastic (DIS). It is a nontrivial task to 
put them together in a Monte Carlo generator of 
events. Most problems arise in the resonance re- 
gion. There are several models of single pion pro- 
duction (SPP) due to resonance production but 
some non-resonant background is also required. 
Then it is necessary to join such a model with 
the DIS part which must be extrapolated far away 
from the kinematical region in which it is reliable. 
In existing MC codes |T| the Rein-Sehgal model 
[2] is most often used to describe the dynamics of 
SPP. It includes contributions from several reso- 
nances of the mass up to 1.8 GeV added in the 
coherent way. The non-resonant background is 
then added incoherently in order to fit the ex- 
perimental data. An important improvements to 
existing MC codes can come from precise JLab 
experimental data from electron scattering exper- 
iments 0. It is known how to construct good 
experimental fits for the structure functions Fi 
and F2 which in the resonance kinematical region 
average over resonance peaks ^ . In the neutrino- 
nucleon cross section the axial structure function 
F3 is also present whose modification cannot be 
deduced from electron experiments. As a first 
guess one can assume that the modification is 
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Figure 1. Fraction of single pion production 
contribution in overall Deep Inelastic Scattering 
cross section 



analogous to that applied to the Fi^2 and investi- 
gate consequences of such assumption. 

This was the starting point for our investiga- 
tion. We shall describe a construction of MC 
code WroNG (WROclaw Neutrino Generator) 
for SPP which incorporates explicit A excita- 
tion model and three exclusive SPP channels ex- 
tracted from the DIS formalism. We focused 
on SPP channels since inclusion of the remain- 
ing dynamics (quasi-elastic channel and more in- 
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Figure 2. Total cross section in the channel 
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Figure 3. Total cross section in the channel 



elastic reactions described by DIS formalism) is 
straightforward. Our code can supplement the 
NUX+FLUKA scheme which does not contain a 
resonance contribution 5 . 

In order to evaluate SPP in the framework 
of the DIS formalism we introduced three func- 
tions (one for each channel) of kinematical vari- 
ables. The functions measure the probability 
that after fragmentation and hadronisation the 
final hadronic state is that of SPP. We obtained 
these functions from the NUX+FLUKA simula- 
tions which are based on the LUND algorithm 6 . 
They turned out to be monotonously decreasing 
functions of the hadronic mass VF, taking values 
in a range from 1 to (see fig. 1). 

The functions we introduced above have been 
used to define the differential cross section for 
SPP in the DIS formalism (there are three iden- 
tical definitions for each reaction channel): 
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where the DIS inclusive cross section. 

W and u) denote hadronic mass and energy trans- 
fer respectively. 



Figure 4. Total cross section in the channel v^p 
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Figure 5. Hadronic mass distribution in the chan- 
/Li~7r+n at = IGeV 
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Figure 6. Hadronic mass distribution in the chan- 
H'TT^p at = IGeV 



nel I'^n 



each reaction channel can be written as: 



In our work we had to solve two problems. 
The first was to join A and DIS contributions. 
The second problem was to add appropriate non- 
resonant contributions. In what follows by DIS 
we mean exclusive SPP channels contained in 
inclusive DIS. We joined two dynamical mech- 
anisms in the cross section expression accord- 
ing to the values of the hadronic invariant mass 
in the kinematically allowed region. The basic 
idea was that for small (i.e. from the threshold 
W = M + m.„ to about 1.4 GeV) values of W the 
dynamics is that of A excitation while for larger 
W the dynamics is that of DIS. In order to make 
the transition smooth we fixed a region from Wi 
to W2 in which the probability to choose either 
of two dynamics changed linearly in the MC way. 

We mimic the "non-resonant" background in 
the region of small values of W by an admixture 
of the DIS contribution. It was done in the MC 
way and amount of the DIS contribution was de- 
scribed by a parameter a a value of which was 
fixed by making a comparison with experimen- 
tal data separately in each exclusive channel. To 
summarize the formula for the cross section in 
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2. RESULTS 



First we present basic ingredients of our con- 
struction: the A excitation model (taken from 
the Marteau model 0) and the model of SPP 
based on DIS. We show the total cross sections 
in three reaction channels (fig. 2-4) with exper- 
imental points taken from papers: [SI- We 
also show distributions of events in the hadronic 
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Figure 7. Hadronic mass distribution in the chan- 
nel Vp^p — > /i~7r+p at El, = IGeV 



Figure 8. WroNG predictions in the channel 
Vf^ji — !■ /i~7r+n and experimental data (total cross 
section) from |H] [HI 



mass for neutrinos of energy 1 GeV (fig. 5-7). In 
the DIS part we did computations for two sets of 
PDF: modifications of GRV94 and GRV98. Dif- 
ferences between them are very small and in what 
follows we restrict ourselves to GRV98 with cor- 
rections. 

In the channel f^n —>■ fi^ir^n the DIS cross 
section is much bigger than for the A excitation 
one (see fig. |2l . If compared with experimental 
data we see that the DIS predictions are above ex- 
perimentally observed while those of the A model 
below. In fig.|31the DIS differential cross section 
of hadronic mass is bigger than of the A in the 
whole kinematical domain. 

In the channel ly^^n /i^7r*^p the situation is 
quite different (see fig. ISJ- At smaller values of 
the neutrino energy E (about 1 GeV) both mod- 
els predict similar values of the total cross section 
close to experimentally measured. For higher val- 
ues of neutrino energy the DIS predictions agree 
with the experimental data while the A excita- 
tion model predictions are too low. In fig. |B1 dif- 
ferential cross sections reveal that while the cross 
sections are similar in the total area below the 
curves, the shapes of the hadronic mass distribu- 
tions of events are very different. In this channel 



we see a nice manifestation of the quark-hadron 
duality: the DIS contribution average over reso- 
nance peak. 

In the channel i^^p /i^vr+p the A excitation 
model predicts much higher values of the cross 
section which are also close to the experimentally 
measured (see fig. 010) • 

We conclude that each channel has its unique 
features and has to be treated independently. We 
also note that quark-hadron duality is seen in only 
one SPP channel. It suggests that modifications 
of PDF we applied in our computations are not 
yet good enough. 

Many choices for Wi , W2 for each channel sep- 
arately were checked in order to find the most 
suitable one. The results did not depend much 
on the choice and we fixed for all three channels: 
W^i = 1.3 GeV and 1^2 = 1-6 GeV. 

Total cross sections depend in a substantial way 
on a. In the channel v^p —>■ fj,~Tr~^p an increase 
of a makes the cross section smaller, in other two 
channels the dependence is opposite. We also 
looked at distributions of events in the hadronic 
mass. We compared our MC results with Kita- 
gaki et al. data [5] because they have the best 



5 




1 1.5 2 2.5 3 1 2 3 4 5 



W(GeV) neutrino energy (GeV) 



Figure 9. WroNG predictions in the 

channelt'pn —>■ fi^n^n and experimental 
data (hadronic mass distribution) from ^ 



Figure 10. WroNG predictions in the channel 
t^^n — /i~7r°p and experimental data (total cross 
section) from |H] [HI 



statistics. We took into account experimentally 
reconstructed spectra of neutrinos and we pro- 
duced samples of events with the same spectrum. 
The number of events has been chosen to be the 
same as in the oryginal experiment so that this 
part of the analysis applied only to the shapes of 
hadronic mass distributions of events. 

In two channels on the neutron the value of a 
determines the height of the A resonance peak. 
The transition region (W^i,W2) can manifest it- 
self as a resonance-like peak at values of W close 
to W2 ■ It is only by chance that such higher reso- 
nance peak can be seen in the experimental data. 
With an increase of a the resonance peak becomes 
lower but at the expense of too many events at 
lower values of hadronic mass and in clear contra- 
diction with experimental data. The best value 
of a is a compromise between two described ten- 
dencies. 

In figs. 8-13 we show a comparison of our best 
qualitative fits with Kitagaki et al. experimental 
data in three reaction channels. We took a = 0.2 
for v^n — > ^"TT+n channel, a — 0.3 for i^^n 
H~Tr'^p channel and a = for Vf^p fi^n^p 
channel. Taking into account the simplicity of 



our construction we find the agreement satisfac- 
tory. 

We plan to improve WroNG in the following 
way: 

• a non-resonant background will be intro- 
duced in theoretically justified way using 
Fogli-Nardulli [E] or Lee-Sato dH results 

• three functions Fspp{W) will be derived di- 
rectly from PYTHIA/JETSET or LEPTO 

• the model will be extended to describe SPP 
on nucleus targets 
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i^„,n —> /x^TT^p and experimental data (hadronic 
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Figure 12. WroNG predictions in the channel 
Vfj^p — > fi~'!T^p and experimental data (total cross 
section) from |H] (H] 
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Figure 13. WroNG predictions in the channel 
Vfj^p — > pT'K'^p and experimental data (hadronic 
mass distribution) from [S] 



